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We present a search for Galactic dark matter (DM) satellites using the Large Area Telescope (LAT). N-body 
simulations based on the ACDM model of cosmology predict a large number of as yet unobserved Galactic DM 
satellites. These satellites could potentially produce 7-rays through the self-annihilation of DM particles. Some 
DM satellites are expected to have hard 7-ray spectra, finite angular extents, and a lack of counterparts at other 
wavelengths. We searched for LAT sources with these characteristics. We found no candidate DM satellites 
matching these criteria in one year of LAT data and interpreted this result in the context of N-body simulations. 



1. Introduction 

We report on a search for dark matter (DM) satel- 
lites via 7-ray emission from weakly interacting mas- 
sive particle (WIMP) annihilation. Specifically, we 
searched for satellites predicted by cosmological N- 
body simulations [3, [2j but lacking counterparts in 
other wavelengths. We selected unassociated, high- 
Galactic-latitude 7-ray sources from both the First 
LAT Source Catalog (1FGL) [3[ and an independent 
list of source candidates created with looser assump- 
tions on the source spectrum. Using the likelihood ra- 
tio test, we distinguished extended sources from point 
sources and WIMP annihilation spectra from conven- 
tional power-law spectra. No candidates were found 
in either the unassociated 1FGL sources or our addi- 
tional list of candidate sources. This null detection 
is combined with the Via Lactea II (VL-II) [l| and 
Aquarius Q simulations to set an upper limit on the 
annihilation cross section for a 100 GeV WIMP anni- 
hilating through the bb channel. 



2. Analysis 

2.1. Data Selection 

Our data sample consisted of 'Diffuse' class events 
with energies between 200 MeV and 300 GeV from the 
first year of LAT data collection (2008 August 8 - 
2009 August 7). We rejected events with zenith angles 
larger than 105° and events taken during time peri- 
ods when the rocking angle of the LAT was greater 
than 47°. This analysis was limited to sources with 
Galactic latitudes greater than 20°, since the Galac- 
tic diffuse emission complicates source detection and 
the analysis of spatial extension at lower Galactic lat- 
itudes. We modeled the diffuse 7-ray emission with 
the standard Galactic (glLiem-v02.fit) and isotropic 



(isotropic_iem-v02.txt) background models. Through- 
out this analysis, we used the LAT ScienceTools ver- 
sion v9rl8pl and the P 6 _V3 -DIFFUSE instrument re- 
sponse functions. 1 

2.2. Source Selection 

The 1FGL contains 1451 high-energy 7-ray sources, 
of which 806 are at high Galactic latitude (\b\ > 20°). 
Of these high-latitude 1FGL sources, 231 are unasso- 
ciated with sources at other wavelengths and consti- 
tute the majority of the sources that were examined 
as potential DM satellites. The 1FGL spectral anal- 
ysis, including the threshold for source acceptance, 
assumed that sources were point-like with power-law 
spectra. This decreased the sensitivity of the 1FGL 
to both spatially extended and non-power-law sources, 
which are characteristics expected for DM satellites. 
To mitigate these biases, we augmented the unassoci- 
ated sources in the 1FGL with an independent search 
of the high-latitude sky [4[ . 

We searched for 7-ray sources using the internal 
LAT Collaboration software package, Sourcelike [5]. 
Sourcelike performs a fully binned likelihood fit in 
two dimensions of space and one dimension of energy. 
For spectral fitting, Sourcelike fits the number of 
counts associated with a source in each energy bin 
independently. The full likelihood is the product of 
the likelihoods in each bin. This calculation has more 
degrees of freedom than that performed by the LAT 
ScienceTool, gtlike, which calculates the likelihood 
from all energy bins simultaneously according to a 
user-supplied spectral model. In this analysis, we used 
11 energy bins logarithmically spaced from 200 MeV 
to 300 GeV. 

Using Sourcelike, we searched for sources in 2496 
regions of interest (ROIs) of dimension 10° x 10° cen- 
tered on HEALPix [6] pixels obtained from an order 4 
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tessellation of the high- latitude sky (\b\ > 20°). Each 
ROI was sub-divided into 0.1° x 0.1° pixels, and for 
each pixel the likelihood of a point source at that loca- 
tion was evaluated by comparing the maximum likeli- 
hood (C) of two hypotheses: (1) that the data were de- 
scribed by the standard LAT diffuse background mod- 
els without any point sources (Ho), and (2) that the 
data were described by the existing model with an 
additional free parameter corresponding to the flux 
of a source at the target location (Hi). Utilizing 
the likelihood ratio test, we defined a test statistic, 
TS = -2ln(C(H )/C(Hi)). 

After creating a map of TS over the high-latitude 
sky, candidate sources with TS > 16 were iteratively 
refit incorporating the normalizations of diffuse back- 
grounds and neighboring point sources as free parame- 
ters. After refitting, source candidates with TS > 24 2 
were accepted for study. Finally, to avoid duplicat- 
ing 1FGL sources, we removed candidate sources with 
68% localization errors overlapping the 95% error el- 
lipse given in the 1FGL. 

Our search of the high-latitude sky revealed 710 
candidate sources, of which 154 are not in the 1FGL 
(36 of these candidate sources were subsequently in- 
cluded in the Second LAT Source Catalog (2FGL) Q). 
We did not expect to recover all of the 806 high- 
latitude 1FGL sources, since the 1FGL is a union of 
four different detection methods and external seeds 
from the BZCAT and WMAP catalogs Q. However, 
since Sourcelike fits each energy bin independently, 
we expected to find source candidates that were not 
included in the 1FGL, either because they had non- 
power-law spectra or they had hard spectra with too 
few photons to pass the 1FGL spectral analysis. We 
sacrificed some sample purity for detection efficiency 
in our candidate source list because stringent cuts on 
spatial extent and spectral shape were later applied. 
We obtained a final list of 385 high-latitude unassoci- 
ated LAT sources and source candidates by combining 
231 unassociated 1FGL sources with 154 non-lFGL 
candidate sources. 

To check for consistency with the source analysis 
of the 1FGL, we performed an unbinned likelihood 
analysis with gtlike assuming that the unassociated 
sources were point-like with power-law spectra. Our 
fitted fluxes and spectral indices are in good agree- 
ment with those in the 1FGL for the 231 unassociated 
1FGL sources. The wide range of fluxes and spectral 
indices spanned by the unassociated LAT sources can 
be seen in Figure [U It is apparent that there are more 
non-lFGL source candidates in this sample with very 
hard spectra (spectral index ~ 1.0) and very low fluxes 
(~ 10 -10 phcm -2 s _1 ). We find that these source can- 
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Figure 1: The distribution of spectral indices and 
integral fluxes from 200 MeV to 300 GeV for the 385 
high- latitude unassociated sources and source candidates. 
The squares are the 231 unassociated sources from the 
1FGL, while the triangles are the 154 additional source 
candidates detected with Sourcelike. The stars are the 
10 representative power-law models in Table HI 



didates are very likely spurious 



2.3. DM Satellite Candidate Selection 

The 7-ray flux from WIMP annihilation in a satel- 
lite can be expressed as 



%IMP 



(1) 



2 Simulations show that 1 in 10 4 background fluctuations are 
detected at TS > 24 when fit with Sourcelike. 0. 



Here, J(^) represents the line-of-sight integral 
through the DM density at an offset angle ip (rela- 
tive to the center of the satellite), while & PP (E) is 
an energy-dependent particle physics factor including 
the velocity- averaged annihilation cross section (<jv) 
(see |8[ for more details). 

While the majority of DM satellites in current simu- 
lations are not spatially resolvable by the LAT, spatial 
extension is an important feature for distinguishing 
large or nearby satellites from point-like astrophysi- 
cal sources. Assuming that the DM distribution of a 
satellite follows an NFW profile [9] with scale radius 
r s at a distance D, the angular extent of the satellite 
can be characterized by the parameter ao = r s /D. 
Approximately 90% of the integrated J-factor comes 
from within the angular radius c^o [lQ[ • 

We used the likelihood ratio test, as implemented 
by Sourcelike, to test sources for spatial extension. 
We defined a test statistic for extension as TS ext = 
TSnfw — TSpoint, where TS po int is the test statistic of 
the candidate source assuming that it has negligible 
extension (ao much smaller than the LAT PSF) and 
TSnfw is the test statistic of the candidate source 
when ao is fit as a free parameter. In both cases, the 
position of the source was optimized during the fit. 

Additionally, we defined a test statistic to select 
sources with 7-ray spectra consistent with WIMP an- 
nihilation into bb (chosen as a representative proxy 
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Table I Values for TS^ t and TS^ ec 



Model Spectral 



Flux (a) 



rpq99 rpq99 



Number Index ( ph cm 2 s 1 ) 

1 0.9 2.0 x 10" 10 

2 0.9 8.0 x 10" 11 

3 1.5 1.1 x 10" 9 

4 1.5 2.0 x 10" 10 

5 2.0 1.2 x 10" 8 

6 2.0 1.2 x 10~ 9 

7 2.5 2.1 x 10" 8 

8 2.5 0.5 x 10" 8 

9 3.0 1.7 x 10~ 8 

10 3.0 1.0 x 10" 8 



6.18 
7.87 
5.09 
14.98 
5.11 
9.63 
6.74 
10.78 
9.81 
11.87 



2.38 
2.46 
4.96 
2.88 
2.24 
4.28 
1.78 
5.66 
2.14 
6.02 



(a) Integral flux from 200 MeV to 300 GeV. 
Note - Cuts for excluding point-like sources and 
sources with power-law spectra. 



for tree- level annihilation spectra). This spectral test 
statistic, TS spec = TS 6 5 — TS pw i, is the difference in 
source TS calculated with an unbinned analysis using 
gtlike assuming a bb spectral model (TS55) and a 
power-law (TS pw i) spectral model. When performing 
our fit, we modeled the candidate source as point-like 3 
and left its flux and the flux of the diffuse backgrounds 
free. In addition to the flux, each spectral model con- 
tained an additional free parameter (the DM mass or 
spectral index). 

We defined cuts to independently eliminate 99% of 
point sources and 99% of power-law sources, labelled 



TS ext and TS ext respectively. To evaluate these cuts 
over the pertinent range of source fluxes and spectral 
indices, we bracketed the unassociated LAT sources 
with 10 representative power-law models (Table U and 
blue stars in Figure [1]). For each of these 10 mod- 
els, we simulated 1000 independent sources at ran- 
dom locations in the high-latitude sky using the LAT 
simulation tool, gtobssim, and the spacecraft point- 
ing history for our one-year data set. To accurately 
incorporate imperfect modeling of background point 
and diffuse sources, we embedded the simulated point 
sources in the LAT data and calculated TS ext and 
TS 



spec* 



We defined TS^ t and TS^ ec for each rep- 



resentative model as the smallest value of TS ext or 
TS spec that was larger than that calculated for 99% 
of simulated power-law point sources. The values of 
TSg Xt and TSgp ec (Table IJ) were calculated indepen- 
dently of each other and the TS > 24 detection cut. 
We used a bilinear interpolation of Table [I] to esti- 
mate the value of TSg Xt and TSgp ec for any point in 
the space spanned by the grid of flux and spectral 



3 Assuming that sources are point-like was found to be a 
conservative way to estimate TS spec 



index. These tests of spatial extension and spectral 
character allow us to select non-point-like and non- 
power-law sources with a contamination of 1 in 10 4 
assuming they are independent. 



3. Results 

We applied the cuts on spatial extension and 
spectral character to select DM satellite candi- 
dates from the 385 unassociated high-latitude LAT 
sources and source candidates. Two of the 
385 unassociated sources, 1FGL J1302.3— 3255 and 
1FGL J2325.8— 4043, passed the cut on spatial exten- 
sion. One of these, 1FGL J1302.3— 3255, also passed 
our spectral test, preferring a bb spectrum to a power- 
law spectrum. However, we do not believe that either 
of these sources is a viable DM satellite candidate for 
reasons discussed below. 

1FGL J1302.3— 3255 was unassociated when the 
1FGL was published but has since been associated 
with a millisecond pulsar by radio follow-up observa- 
tion [n|. The other source, 1FGL J2325. 8-4043, has 
a high probability of association with two AGN [121 ], 
though it did not meet the association criteria of 
the 1FGL. Cross checking against the 2FGL, two 
sources were found within 0.5° of the location of 
1FGLJ2325. 8-4043 0. In one year of data, these 
two sources could not be spatially resolved, but their 
existence was enough to favor an extended source hy- 
pothesis. 

Since 1FGL J1302.3— 3255 was associated with a 
pulsar and 1FGL J2325.8— 4043 did not appear to 
be truly extended, we conclude that there were no 
unassociated, high-latitude spatially extended 7-ray 
sources in the first year of LAT data. Thus, accord- 
ing to the criteria defined in Section [2j no viable DM 
satellite candidates were found. 



4. Discussion 

Using the detection efficiency of our selection, the 
absence of DM satellite candidates can be combined 
with the Aquarius and VL-II simulations to constrain 
a conventional 100 GeV WIMP annihilating through 
the bb channel. We calculated the probability of de- 
tecting no satellites from the individual detection effi- 
ciency of each satellite in the realization. By increas- 
ing the satellite flux until the probability of detecting 
no satellites drops below 5%, we set a 95% confidence 
upper limit on (<jv). 

4.1. Detection Efficiency 

The detection efficiency of our selection is defined 
as the fraction of true DM satellites that would pass 
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Table II Detection efficiency 



Flux (a) 




Extension 




( phcm -2 s _1 ) 


0.5° 


1.0° 


2.0° 


0.2 x 10" 8 


< 0.05 


< 0.05 


< 0.05 


0.5 x 10" 8 


0.16 


0.28 


0.31 


1.0 x 10~ 8 


0.74 


0.76 


0.83 


2.0 x 10" 8 


0.99 


1.0 


0.99 


5.0 x 10" 8 


1.0 


1.0 


1.0 


(a) Integral flux from 200 MeV to 300 GeV. 
Note - Satellite detection efficiency for a 
100 GeV WIMP annihilating to bb. 



the cuts in Section [2] and was calculated from Monte 
Carlo simulations. The efficiency for detecting a DM 
satellite depends on spectral shape (i.e., DM mass and 
annihilation channel), flux, and spatial extension. For 
a 100 GeV WIMP annihilating through the bb channel, 
we examined the efficiency for satellites with fluxes 
in the range of the unassociated high-latitude LAT 
sources and angular extents in the range to which the 
LAT is currently sensitive (Table HI]). 

For each set of characteristics listed in Table [III 
we simulated 200 DM sources with NFW profiles 
and bb spectra from a 100 GeV WIMP. These simu- 
lations were embedded in LAT data at random high- 
latitude locations, and Sourcelike was used to com- 
pute TS ext , TS spec , and the detection TS for each. 
The satellite detection efficiency was computed as 
the fraction of satellites with Sourcelike TS > 24, 
TS ext > TS g Ji t , and TS spec > TS^ ec . The first re- 
quirement was included as a proxy for the efficiency 
of the source finding algorithm. To expedite the gen- 
eration of this table, we found the flux value with 
efficiency < 0.05 and conservatively set the efficiency 
for sources with less flux to 0. 

4.2. Upper Limits 

Realizations of the Galactic DM satellite popula- 
tion were created by selecting 6 maximally separated 
vantage points 8.5 kpc from the center of the VL- 
II simulation and each of the 6 Aquarius simulations. 
For each of the 6 x (1 + 6) = 42 "visualizations" of VL- 
II and Aquarius, we calculated the 7-ray fluxes of all 
satellites for a given (av) using Eqn. (pQ). With these 
fluxes and the true spatial extension for each satellite, 
we performed a bilinear interpolation on Table [IT] to 
determine the detection efficiency for each satellite. 
The probability that the LAT would observe none of 
the satellites in visualization i is 

P i ((av))=U(l-e hJ ((vv))) (2) 

3 

where e^j is the detection efficiency for satellite j in 
visualization i. Because there is no reason to favor 



any one visualization, we calculated the average null 
detection probability over the N = 42 visualizations 

as 

1 N 

^(M) = ^E P *(M) (3) 

i 

To set an upper limit on the DM annihilation cross 
section, we increased (av) until the probability of a 
null observation was < 5%, i.e. P < 0.05. This corre- 
sponds to 95% probability that, for this (av), at least 
one satellite would have passed our selection criteria. 
Using this methodology, the LAT null detection con- 
strains (av) to be less than 1.95 x 10 _24 cm 3 s _1 for a 
100 GeV WIMP annihilating through the bb channel. 
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